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Abstract 

There has been much recent research into polarizing an antiproton beam, in- 
stigated by the recent proposal from the PAX (Polarized Antiproton experiment) 
project at GSI Darmstadt. It plans to polarize an antiproton beam by repeated 
interaction with a polarized internal target in a storage ring. The method of po- 
larization by spin filtering requires many of the beam particles to remain within 
the ring after scattering off the polarized internal target via electromagnetic and 
hadronic interactions. We present and solve sets of differential equations which de- 
scribe the buildup of polarization by spin filtering in many different scenarios of 
interest to projects planning to produce high intensity polarized beams. These sce- 
narios are: 1) spin filtering of a fully stored beam, 2) spin filtering while the beam 
is being accumulated, i.e. unpolarized particles are continuously being fed into the 
beam, 3) the particle input rate is equal to the rate at which particles are being 
lost due to scattering beyond ring acceptance angle, the beam intensity remaining 
constant, 4) increasing the initial polarization of a stored beam by spin filtering, 5) 
the input of particles into the beam is stopped after a certain amount of time, but 
spin filtering continues. The rate of depolarization of a stored polarized beam on 
passing through an electron cooler is also shown to be negligible. 

PACS: 13.88.+e Polarization in interactions and scattering - 24.70.+S Polarization phe- 
nomena in reactions - 25.43.-t-t Antiproton-induced reactions - 29.27.Hj Polarized beams 

1 Introduction 



The availability of a stored polarized antiproton beam offers many new opportunities 
for studying the structure of hadronic matter. The collision of polarized protons with 
polarized antiprotons at relatively high luminosity, for example, enables an evaluation 
of the transverse polarization of a quark within a proton. Such a determination would 
provide an important element of a QCD description of the partonic structure of a nucleon. 
Transversity can only readily be accessed through a Drell-Yan process induced by the 
scattering of appropriately polarized protons and polarized antiprotons [H [21 [3l HJ [5] . 

Another area of study that would be opened up as a consequence of having an intense 
beam of polarized antiprotons relates to the time-like electromagnetic form factors of a 
proton, particularly their phases, since the form factors take complex values in the time- 
like region above threshold j6j[7]. The expectations of perturbative QCD could be probed 
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in some detail with such form factor measurements, in addition to further results arising 
from the hard scattering of polarized antiprotons and polarized protons at high Q 2 [5]. 

The buildup of the polarization of an antiproton beam by spin filtering off a polar- 
ized internal target (PIT) has been described by differential equations [HJ [9j [TOj [HI [12] 
involving the spin observables for antiproton-electron and antiproton-proton elastic scat- 
tering, building on the earlier work of refs. [T3], [HI [15] . This is of interest to the PAX 
collaboration who plan to polarize an antiproton beam by repeated interaction with a 
polarized internal target in a storage ring [5], [16] . Relativistic expressions for spin observ- 
ables due to single photon exchange in elastic spin 1/2 - spin 1/2 particle collisions have 
been presented recently [IT]. We use these spin observables to investigate the buildup of 
polarization of a beam of spin 1/2 particles by spin filtering. 

The paper is organized as follows: the method of polarization buildup by spin filtering 
is described in sect. [5] and later sections each investigate a particular scenario of spin 
filtering. Firstly spin filtering of a fully stored beam is described in sect. Eland then spin 
filtering as the beam is being accumulated is described in sect. [H In sect. [5] a system with 
constant beam intensity is analyzed, i.e. the input rate exactly balances the rate of loss 
of beam particles due to scattering out of the ring. Polarization buildup of a stored beam 
which is initially polarized is the subject of sect. EJ and finally a system where particles 
are fed into the beam for a limited time after which spin filtering continues is treated in 
sect. [7J The paper is summarized in sect. 

2 Spin filtering 

Spin filtering is the idea that a beam in a storage ring can be polarized by repeated 
interaction with a polarized internal target because of a difference between the cross- 
section for scattering of particles with their spins aligned parallel and the cross-section 
for scattering of particles with their spins aligned anti-parallel. It was first proposed by 
P. L. Csonka in 1968 [18J. A ring acceptance angle 9^ divides the range of scattering angle 
9 into two distinct regions, a particle scattered an an angle below the acceptance angle 
remains "in" in beam and a particle scattered at an angle greater than the acceptance 
angle gets scattered "out" of the beam. There is also a minimum scattering angle # m in> 
corresponding to the average transverse electron separation for a pure electron target or to 
the Bohr radius of the atoms in an atomic target, below which scattering is prevented by 
Coulomb screening. The two physical processes that contribute to polarization buildup in 
spin filtering are: (a) spin selective scattering out of the ring, and (b) selective spin-flip. 
Thus particles in one spin state may be scattered out of the beam, or have their spin 
flipped while remaining in the beam, at a higher rate than particles in the other spin 
state. Thus over time one spin state is depleted more than the other leading to a beam 
polarization. 

As the beam polarization increases the beam intensity decreases because particles are 
being continuously lost from the beam. Thus there is a trade off between beam polariza- 
tion and beam intensity that is characteristic of spin filtering. The Figure Of Merit (FOM), 
defined as the polarization squared times the beam intensity FOM(r) = P 2 (r) iV(r), 
where r is the time elapsed since spin filtering began, provides a measure of the quality of 
the polarized beam taking this trade off into account. We are investigating continuously 
inputing unpolarized particles into the beam, during spin filtering, to compensate the loss 
of beam intensity due to scattering out of the ring. 
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Note some treatments of spin filtering investigate a scenario where no particles are 
scattered out of the beam, i.e. the maximum scattering angle for the process is less than 
the ring acceptance angle, which is the case for antiprotons scattering off electrons in an 
atomic target [8j [9j [TOj [HI [15] and for antiprotons scattering off a co-moving beam of 
electrons or positrons [12]. In these scenarios only selective spin- flip can contribute to 
polarization buildup, and one avoids the problem of decreasing beam intensity. The low 
density of the targets currently available causes the rate of polarization buildup using 
these methods to be slow, but the enhanced cross-sections at low energies suggested in 
refs. pL2~, 19] may compensate this difficulty A comparison of all recent treatments of 
spin filtering, including a derivation of the polarization evolution equations treated here, 
is presented in ref. |20j. 

Spin filtering was demonstrated experimentally for polarized protons in 1993 by the 
FILTEX experiment at the Test Storage Ring in Heidelberg [2TJ . The spin filtering method 
is at the heart of the recent PAX proposal to generate a polarized beam of antiprotons in 
the HESR ring of FAIR at GSI Darmstadt [16]. 



3 Polarization buildup of a stored beam 



In this and the next few sections we describe systems of differential equations that model 
the buildup of polarization of a beam by spin filtering in a storage ring. First a system 
where the beam has already been accumulated in the storage ring is analyzed. The 
number of particles in the beam initially is Nq and in this case the number of particles in 
the beam decreases continuously because of scattering out of the ring. When circulating 
at frequency u, for a time r, in a ring with a polarized internal target of areal density n 
and polarization Vt oriented normal to the ring plane, (or longitudinally with rotators) 
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describes the rate of change of the number of beam particles N(t) = N^(t) + iV|(r) and 
their total spin J(r) = iVf(r) — ATj_(r) [H[9]. These coupled differential equations involve 
angular integration of the spin observables presented in ref. [17] over the following ranges, 
as seen in table 1. The "in" subscript refers to particles that are scattered at small 
angles < 9 acc remaining in the beam, and the "out" subscript refers to particles that are 
scattered out of the beam. Thus the integrals over scattering angle 9 are labeled "in" 
where the range of integration is # m i n < 9 < 9 acc , "out" where the range of integration is 
$acc < 9 <n and "all" = "in" + "out" where the range of integration is 9 m \ n < 9 < 7r; as 
seen in table 1. I — da / dQ is the spin averaged differential cross-section and A, K and D 
are the double spin asymmetry, polarization transfer and depolarization spin observables 
respectively as calculated in ref. [17]. All cross-sections and spin observables contributing 
to spin filtering are azimuthally averaged, due to the geometry of the scattering, where the 
scattering plane can be at any azimuthal angle. Hence single spin observables, for example 
the analyzing power, do not contribute to the polarization evolution equations because 
they vanish when azimuthally averaged. The parameters in the matrix of coefficients of 
eq. (13. ip depend on the state of the target polarization, i. e. longitudinal or transverse, as 
seen in table 1. 

We now solve the set of polarization evolution equations presented in eq. (13. ip . The 
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eigenvalues of the matrix of coefficients are found to be 



\ 1 = -nv (J out + L in + L d ) , 

A 2 = -nv (J out + L in - L d ) , (3.2) 
where the discriminant Ld of the quadratic equation for the eigenvalues is 



V2A out (A M - K- m ) + Lf n , (3.3) 

and L in = (J in — D in ) / 2 is a loss of polarization quantity. Note that J out , L in and L d 
are all non-negative. As a consequence the eigenvalues are non-positive and Ai < A 2 < 0. 
When there is no scattering out of the ring all of the "out" integrations are zero, and 
one finds Ai = —2nvL d and A 2 = 0. 

Now enforcing the initial conditions N(0) = N the total number of particles in the 
beam initially, and J(0) = =>■ iVf(O) = iV|(0) = Nq / 2 i.e. initially the beam is unpolar- 
ized, gives the solutions: 

AT , s [e XlT (L d -L in ) + e A ^ (L d + L in )] N 

N (V = ~ 77J , (3-4) 

ZL d 

J(V = ttj • (3-5) 

Z L d 

The time (r) dependence of the polarization of the beam is given by 

V ( T ) = = -(A aU -K- m )V T 

N(t) L in + Ld coth (La n v r) 

The expression for Vir) is proportional to Vt which confirms that if the target polarization 
is zero there will be no polarization buildup in the beam, as was required by eqs. (13.11) . 
The approximate rate of change of polarization for sufficiently short times, and the limit 
of the polarization for large times are respectively: 



dV 
dr 



nvV T {A all - K in ) , (3.7) 



P max = lim V(t) =-V t A ; U - f in . (3.8) 

r^oo L in + L d 

For pure electromagnetic scattering the double spin asymmetries equal the polarization 
transfer spin observables [17] . thus one can simplify the above equations using A- in = K in , 
A on t = K out and A aU = K aU ; hence A a n - K in = K out . 

Eqs. (13. ip neglect a spin tune effect |9j, induced by a pseudomagnetic field generated 
by the target polarization dependent real part of the forward scattering amplitude. The 
effect which is expected to be comparatively small comes from the difference between the 
spin elastic helicity amplitudes, which are known for electromagnetic scattering from QED 
[T71 12"3"] . but not well known for hadronic scattering even though they could be evaluated 
from dispersion relations for the appropriate spin dependent total cross-sections. 
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Transverse polarization requires 



Longitudinal polarization requires 
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Table 1: The entries in the system of equations for polarization buildup involve angular 
integration over the spin observables presented in refs. [ITJ [TJJJ [22] • X, Y and Z are the 
coordinate axis where the beam is moving in the positive Z direction. The minimum value 
for 9 (0 m i n ) relates to the average transverse electron separation for a pure electron target 
and to the Bohr radius for an atomic gas target, and 9 acc is the ring acceptance angle. 



3.1 Beam lifetime and figure of merit 



The beam lifetime r*, the time taken for the beam intensity to decrease by a factor of e, 
can be obtained from eq. (13. 4p . One finds 



1 
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(3.9) 



where J out + L in ^> L d when there is scattering out of the ring. 

The Figure Of Merit (FOM) provides a measure of the quality of the polarized beam, 
and is given by 

FOM(r) = V 2 {r) N{r) J " (r) 



N{t) 



(3.10) 



The figure of merit for the above case is 
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Maximizing the figure of merit gives the optimum polarization buildup time, taking into 
account the trade-off between decreasing beam intensity and increasing beam polarization. 
Solving d FOM / d r = yields T op ti mum ~ 2 / n v ( J out + L m ), approximately twice the 
beam lifetime. Thus the optimum time for polarization buildup is twice the lifetime of 
the beam, as mentioned in ref. [T6l. 
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4 Accumulation of antiprotons in the ring 



In the discussion so far we have only considered polarizing an antiproton beam when 
the beam is already accumulated in the storage ring. The PAX collaboration plans to 
obtain their antiproton beam by collecting the produced antiprotons from high energy 
interactions of protons on targets of light nuclei, such as Beryllium. The antiprotons 
will be continuously fed into the storage ring at a fixed rate and accumulated, hence 
increasing the beam intensity, allowing for a greater luminosity in an experiment. The 
PAX collaboration estimates the production rate of antiprotons as being 10 7 per second 
[5]. Since 10 11 antiprotons are required in the storage ring, antiprotons will be fed into 
the storage ring at a rate of 10 7 per second for 10 4 seconds [5]. 

We now consider a system where spin filtering occurs as the antiprotons are being fed 
into the ring. The original system of equations must be amended to account for this con- 
stant accumulation. The effect will be to add a term (3 to the diV(r) /dr equation, where 
(3 is the constant rate at which antiprotons are fed into the ring; while the d J(r) / d r 
equation remains unchanged. The initial conditions are iV(0) = iVo , which will be set to 
zero in sect. 14.11 an d J{0) = 0. The new system of differential equations is 

diV(r) 



dr 

dJ(r) 
dr 



nu [J out N(t) + V T A ont J(r)\ + (3, 
nu [V T (A aU -K in ) N(t) + (Jan-Din) J(r) 



(4.1) 
(4.2) 



By differentiating eq. f)4.2p with respect to r and substituting in eq. (14.11) one obtains 
an inhomogeneous second order linear differential equation with constant coefficients for 
J(r): 



d 2 JO 



(Ai + A 2 ) 



dJ(r) 



dr 2 

the solution of which is 

2 L d Xi A 2 



dr 



AiA 2 J(r) = -nvV T (A all -K in ) (3. 
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Differentiating eq. (14.41) with respect to r and substituting into eq. (14.21) gives an expression 
for N(t) : 

N(t) = \ [A 2 (X 1 N + /3) (L d -L- m ) e AlT (4.5) 



+ A X (X 2 N + (3) (L iQ + L d ) e X2T + (3 (J all - An) (A 2 -A x )] . 

As a consistency check one can see that these solutions for J(r) and N{r) satisfy the 
initial conditions J(0) = and A^(0) = Nq, and in the particular case when (3 = the 
above expressions reduce to the solution of the homogeneous system presented in sect. [3j 
Dividing J(r) by N(t) we obtain an expression for the polarization V(t) as a function of 
time, 
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When the particle input rate is zero (i.e. (3 = 0) the above equation simplifies to 
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(4.7) 



which is the solution of the homogeneous case presented in eq. (13.61) . 

Using a Taylor Series expansion we find that the approximate initial rate of polarization 
buildup for each of these cases (N Q ^ with (3 ^ and iVo = with (3 ^ 0) is the same 
as in the homogeneous case (N ^ with (3 = 0): 

<±Z ~ -nuV T (Aau - K in ) . (4.8) 
dr 

The maximum polarization achievable is the limit as time approaches infinity: 

-V T (A all -K hl ) 
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lim V(t) 



/all" Di 



(4.9) 



which is independent of both Nq and (3, however note that in taking this limit we used 
the fact that (3^0. If (3 was equal to zero then the maximum polarization achievable 
would equal that from the homogeneous case; as can be easily seen from eq. (14.61) remem- 
bering that Ai < A 2 < 0. Thus for the complete case there are just two values of the 
maximum polarization, one for (3 = and one for all (3^0. The figure of merit for this 
inhomogeneous case is: 
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(4.10) 



where for convenience we have defined the two constants Ci = X 2 (Xi N + (3) and c 2 = 
Ai ( A 2 iVo + (3). Note the FOM will not have a maximum in finite time if the accumulation 
rate (3 is high enough to make the beam intensity a constant or increase with time. If this 
happens the FOM will increase monotonically. 



4.1 No initial beam 

Of interest is the particular case when iVo = 0, i.e. there are no particles in the beam 
initially. In this case the above solutions simplify to 

(3V T ( Aga — K n 



J(r) 
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V(t) 



2 L d Ai A 2 
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Interestingly the f3 dependence of V(t) vanishes in this case, i.e. the polarization buildup 
rate is independent of the rate at which antiprotons are fed into the ring, if there are no 
particles in the beam initially. But we have used the fact that (3 7^ to obtain the above 
result. We should note the obvious physical fact that if N = and /3 — i.e. there are no 
particles in the beam initially and no particles are fed into the beam, then there will never 
be any particles in the beam; so measuring the beam polarization is meaningless. The 
accumulation rate f3 will effect the figure of merit as it greatly effects the beam intensity 
N(t). 

We can summarize the results for the polarization buildup as 

'0 for 13 = & N = 
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5 Constant beam intensity 

In this case the accumulation rate is set specifically so that extra particles are fed into the 
beam at such a rate so that the beam intensity is kept constant, i. e. fed in at such a rate to 
balance the rate at which particles are scattered out of the beam. The system of equations 
is much simpler in this case. Here N(t) = N is a constant, hence dA^(r)/dr = 0, and 
the J(t) equation becomes a 1st order linear ODE with constant coefficients 

dJ(r) 



dr 



+ nv (J aU - D in ) J(t) 



nv (A all -K- m ) V T Ni 
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(5.1) 



and imposing the initial conditions iV(0) = N and J(0) = one obtains the solution 
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Now the polarization as a function of time can be presented 
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To find the maximum polarization achievable, i.e. the limit as time tends to infinity, we 
note that I a i\ > D in thus —nv{ 7 a n — D in ) < and hence one obtains 

-V T (A ati -K in ) 



= lim V(t) = 



I all — Dh 



(5.4) 
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which is the same as in the inhomogeneous case when /3^0. 

The initial rate of polarization buildup can be obtained by expanding V(t) as a Taylor 
expansion in nvr. Assuming nvr is small we neglect terms of second or higher order 
giving 

~ _ nv (A all -K in ) V T , (5.5) 

a r 

as it was in the homogeneous case presented in eq. (13. 7p . The figure of merit in this case 
is easily obtained 

FOM(r) = NoV r Uaii-^in) 2 [ 1 _ c -»*(/.n-i> ta )T]*_ (5 . 6) 

{I all — Din) 

and increases monotonically with time. 



5.1 Approximating the critical input rate 

The accumulation rate needed to keep the beam intensity constant is important, as this 
critical rate divides the solution of the system into two physically distinct cases. Smaller 
accumulation rates than this critical value cause the beam intensity to decrease, hence 
the FOM will have a maximum in finite time. Larger values than the critical value cause 
the beam intensity to increase continuously, hence the FOM will increase monotonically. 
We can see from eq. H3.4[) that N(t) does not decrease linearly with time r. So the 
accumulation rate needed to keep the beam intensity constant, say /(r), will not be 
linear in r. We now derive the function /(r) and obtain a linear approximation to 
it, which can be used in the inhomogeneous case treated in sect. HI We must solve 
N(r) = N hom (r) + f(r) = N , i.e. 

ZL d 

which leads to 

f( T ) = [2L d -e AlT (L d -L- m )-e x * T (L d + L- m )] . (5.8) 

A linear approximation /l(t) to f(r) can be found by Taylor expanding the exponentials 
to order r (which is valid since nvr is small), to obtain 

f L (r) = nvN Q I out r, (5.9) 

which is in the linear form + f3 c r where f3 c = n v No I ont is the critical value of (3 which 
when added to the diV(r) / dr differential equation in sect. H] approximately makes the 
beam intensity constant. 



6 Stored beam with initial polarization 

We now solve the homogeneous system where the initial polarization is not zero. This will 
be used if two methods of polarizing antiprotons are combined, i.e. if antiprotons were 
produced with a small polarization by some other method and one wanted to increase 
that polarization by spin filtering in a storage ring, where the luminosity could also be 



9 



increased. In this section the beam has been stored and there is no further input of 
particles into the beam. 

The system of differential equations, eigenvalues and eigenvectors are the same as 
sect. |3l but one of the initial conditions is different. The new initial conditions are 
N(0) = N > the total number of particles in the beam initially, and J(0) = Jo 7^ 
=> iVf(O) 7^ iV|(0) i.e. initially the beam is polarized. Note that since the number of 
particles in the beam with one particular spin state must not be greater than the total 
number of particles in the beam the bound | Jq| < is respected. A negative value for 
Jo simply implies that the antiproton beam is initially polarized in the opposite direction 
to the polarization direction of the target. Enforcing these initial conditions leads to the 
solutions 

(J V T A ont - N L in ) (e x ^ - e x ^) + N L d (e x ^ + e x ^) 
N{t) = — (6.1) 

T( > [jV P T (^aii-^in) + JpL in ] (e x ^-e x ^)+J L d (e x ^ + e x ^) 

ZL d 

which reduce to the solutions of the original homogeneous system presented in eqs. (13. 4| 13.5ft 
when Jo — > 0. The beam lifetime is the same to leading approximation as in the homoge- 
neous case when J(0) = 0. Dividing J(r) by N(t) provides the time dependence of the 
polarization of the beam 

v( x L d J - tanh (L d nur) [L in J + N Q V T ( A aU - K in )] 
L d N + tanh ( L d n v r ) [ L in N Q - J V T A out ] 

Denoting the initial polarization V(0) = Jq / Nq = Vq the above can be written as 

v( , L d V - tanh (L d nur) [ L in V Q + V T ( A & \\ - K iQ )] 
[T) L d + ttmh(L d nvr) [L in -V V T A out ] ' l ' j 

The approximate rate of change of polarization for sufficiently short times is found by 
Taylor expanding to first order in r 

^ w nu { [A out V 2 -(A an -K in )] V T -2V L in } . (6.5) 

The limit as time goes to infinity of V(r) in eq. (16.41) is 

Inn V(r) = MLL^I^lIM^I , (6 . 6) 

r ^°° ( Jin + L d ) - A ou tVQ V T 

which of course agrees with the earlier maximum polarization if J = {i.e. Vo = 0). The 
figure of merit for this case is: 

FOM(r) = V 2 (r)N(r) = ^ = (6.7) 

N {[V T (A &ll -K in )+V L in } ( e x ^ -e x ^) +V L d (e x ^ + e x ^)} 2 
2L d {(V V T A out - Lin) (e x ^ - e x ^) + L d (e A ^ + e A ^)} 
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6.1 An unpolarized target 

A special case of this system deserves comment. Given that the beam is initially po- 
larized what happens if the target is unpolarized? One would imagine that the beam 
polarization should decrease and eventually reach zero. We now analyze the equations 
of sect. [6] when the target is unpolarized (i.e. Vt = 0) and use the fact that L d = 
a/ V£ A out (A a ii — K m ) + Lf n = L in when Vt = to obtain the beam polarization as a 
function of time 

V(t) = P e (Al ~ A2)r , (6.8) 

which is an exponentially decreasing function of r for Ai — A2 < 0. The beam polarization 
will not decrease in the special case of Ai — A 2 = 0, but this only happens when L in = 0, 
i.e. when there is no depolarization. The special case of Ai — A 2 = 0, which does not 
lead to polarization buildup as seen from eq. ( 16. 21) . would be avoided by any experimental 
effort, thus is omitted from the rest of the discussion. 

The limit of beam polarization for large times when Vt = is 

lim V{t) = 0. (6.9) 



» 00 



Thus, as expected, if the beam is initially polarized and the target unpolarized then 
the beam polarization will decrease with time and eventually the beam polarization will 
reduce to zero. Thus a beam cannot gain polarization from an unpolarized target by spin 
filtering. 

The figure of merit in this case simplifies to FOM(r) = NoVq e^ 2Xl ^ X2 ^ T which is a 
monotonically decreasing function of r. One can derive a polarization half-life in this case, 
the time taken for the polarization to decrease by a factor of 2, by solving V{r) =Vq/2 
to obtain 

In 2 In 2 . 
n = = — . 6.10 

3 A 2 -Ai 2nuL d v ; 

This scenario occurs in an electron cooler, a device used to focus the beam in many storage 
rings. The beam passes through a co-moving beam of unpolarized electrons with low 
transverse momentum, in order to dampen the transverse momentum of the antiprotons 
in the stored beam. But the low electron areal densities in cooler beams, where typically 
n m 10 -18 fm~ = 1CT 19 mb~ , causes the polarization half-life to be very large. Thus 
this depolarization effect is negligible. 

It has recently been suggested that the positron-antiproton polarization transfer ob- 
servable is very much enhanced at low energies [19] . This enhancement is the basis of the 
recent proposal by Walcher et al. to polarize an antiproton beam by repeated interaction 
with a co-moving polarized positron beam in a storage ring [12J. All of the antiprotons 
remain within the beam in this scenario and one avoids the problem of the antiprotons 
annihilating with protons in an atomic gas target. This large enhancement of the polar- 
ization transfer observable at low energies is due to the unlike charges of the positron and 
antiproton, and does not occur for the like charges case of antiproton-electron scattering. 
Thus this does not effect the conclusion that depolarization of an antiproton beam in an 
electron cooler is negligible. 

6.2 A critical value for the target polarization 

The beam polarization will also decrease for low values of the target polarization. In fact 
there is a critical value of the target polarization Vt which keeps the beam polarization 
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constant. If the target polarization is above this critical value the polarization of the 
beam will increase, and if the target polarization is below this critical value the beam 
polarization will decrease. The critical value is obtained by solving 

rp,s L d V - tanh(L d n?/T) [L in V + V T (A aU - K in )] 
[T) ' L A + tanh(L d nz/r) [L in - V V T A out ] °' l ' J 

for Vti where the time dependence will cancel leading to 

critical _ 2 Pq L j n 

~ VZA mt -(A all -K in )- [b - iZ) 

For target polarizations below this critical value the maximum beam polarization occurs 
at time r = 0, and for target polarizations above this critical value the maximum beam 
polarization occurs at large times r — > oo. 



7 Particles fed in for a limited time 

The Heaviside step function could be used in the system of equations to explain the case 
of particles input into a beam for a certain amount of time after which the input is turned 
off and no more particles are fed into the beam, but spin filtering continues. This scenario 
is under consideration by the PAX collaboration [5]. The Heaviside function is a piecewise 
continuous function which is zero in one region and one everywhere else, it is used in many 
mathematical modeling problems to describe an external effect turned on or off after a 
certain duration of time. In our case this external effect is the input of particles into the 
beam. The Heaviside function is defined as 



H(t-t, 



c 



if T < T c 
if T > Tc 



and is used to describe an external effect turned on at time r c , but we require an external 
effect on initially and turned off at time r c , thus we need 



H(t)-H(t-t c 



if < T < T c 

if r > T r 



I c 



Note that in the special case when r c = 0, [H (r) — H (r) ] = 0, as this is in the second 
region. This describes a physical situation where particles are being fed in for zero seconds, 
which is the same as saying no particles are fed in, so the second order ODE should be 
homogeneous in this case, which it is. The extra term to add to the d N / dr equation 
to account for particles being fed in at a constant rate (3 per second for r c seconds after 
which the input is switched off is (3 [H (r) — H ( r — r c ) ] . 

The equations are now broken into two pieces, i.e. discontinuous, but piecewise con- 
tinuous. The solutions will be broken into two regions < r < r c and r > r c , where the 
solutions in the region < r < r c should equal those in the inhomogeneous case presented 
in sect. SJ The initial conditions will be N(0) = and J(0) = 0, thus J'(0) = 0. The 
Heaviside function above is included in the second order ODE for J(r) to obtain: 

& 2 J{r) f , | \ \ d J(r) 
dr dr 

= -nvV T (A^-K-n) p [H(t)-H(t-t c )] , (7.1) 



12 



which can be solved by Laplace Transform methods giving 
j (t) f A i (l-e A -) + A 2 (e^-1) 



Ci 



if < r < r„ 



if r > T n 



(7.2) 



(7.3) 



Ai ( e A2(r - Tc) - e A2T ) + A 2 ( e AlT - e Al ( T ~ Tc )) 

where for convenience we have defined the constant factor 

c _ (3V T (A aU -K in ) 
2 L d Ai A 2 

One sees from eq. (17.21) that J(r) = for all r when r c = 0, which is physically reasonable 
as there are never any particles in the beam if r c = 0. Also the Vt factor in C\ indicates 
that J(r) will always be zero if Vt = {i.e. if the target is unpolarized). It can also be 
seen that the complete solution in the region r > r c is the combination of the solution in 
the region < r < r c and an additional part dependent on r c ; which is 



Ci [A x ( 



, A 2 ( T-T C ) 



- 1) + A 2 (1 - 



,Ai (t-t c ) 



)] 



and immediately one sees that when r = r c this additional part vanishes. So when r = r c 
i.e. at the boundary between the two regions, the two solutions match. Thus the solution 
for J(r) is continuous as expected. The expression for J(r) in the first region < r < t c 
of eq. (17.21) is equal to the solution of the inhomogeneous system presented in eq. (14. lip . 
A similar analysis as that done for J(r) reveals the second order ODE for N(r) 

d 2 N{r) dN{r) 



dr 2 



Ai + A 2 . 



n v 



'all 



dr 
-D h 



+ X 1 X 2 N{t) 
) /3[H(t)-H( t 



)-H{t-t c )\. 



(7.4) 



The initial conditions are N(0) = and N'(0) = (3 [ 1 — H ( — r c ) ], the latter of which 
deserves comment. The rate N'(0) should be f3 when r c > and when r c = 0, cor- 
responding to the case of particles being fed in for zero seconds, i.e. no particles fed in. 
Note that H ( — r c ) = 1 when r c = and H ( — r c ) = when r c > 0, and note physically 
that r c , the duration for which particles are fed into the beam, cannot be negative. The 
Heaviside function in this initial condition forces the solution for N(t) to be split into 
three regions. On solving by Laplace Transform methods one obtains 



N(t) 

Co 



Ai 



1 



Lin + La 



'all 



A, 



, A2 T 



+ A 2 



Ai 



, A 2 ( T-T C ) _ 



+ A 2 



L m 


-L d 




-D hl 


L m 


+ L d 


Ial\ 


-D hl 


L m 


-L d 


-fall 


-D m 



, Ai r 



if < T < T r 



g Al T _ g Al ( T-T C ) 



(7.5) 



if T > T c > 



if r c = 
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where again for convenience we have defined a constant factor 



C 2 = — - - - — - (7 6) 

2LdAiA2 

Again one sees that the complete solution in the region r > r c > is the combination of 
the solution in the region < r < r c plus an additional part dependent on r c ; which is 

C 2 [Ai ( e A2(r - rc) -l) + A 2 (l- e Al(r - Tc) )] . 

Immediately we see that when r = r c this additional part vanishes, thus the solution for 
N(t) is continuous. The solution in the first region < r < r c is equal to the solution 
from our inhomogeneous case presented in eq. (14.121) . and it satisfies the initial condition 
N(0) = 0. 

We now present results for the polarization V(r) = J(r) / N(r) as a function of time, 
in both regions. The polarization is undefined when there are no particles in the beam, 
thus we need not treat the case r c = 0. As expected in the < r < r c region V(t) 
equals the solution of our inhomogeneous case presented in eq. (14.131) . and in the region 
r > t c > one finds 

V(t) = (7.7) 



V T Maii-^in) [Aie A2T (e~ x ^ - l) +A 2 e 



Ai t ( \ _ g— Al r c 



Aie A ^[(/ all _ J D in ) e -A 2 r c _( L . n + Ld) ] + A2e A 1 r [ (L . n _ Ld ) _ ^ IsLll _ D . n ) e -x 1Tc y 

The approximate initial rate of polarization buildup and the maximum polarization 
achievable will both reside in the < r < r c region, and thus will be identical to those 
presented in sect. HJ This is because the maximum polarization achievable occurs when 
the input rate is never switched off, i.e. in the < r < t c region. 



8 Summary 

We have presented and solved systems of differential equations that describe the buildup 
of polarization of a beam by spin filtering in a storage ring, in various different scenarios. 
These scenarios are: 1) spin filtering of a fully stored beam, 2) spin filtering while the 
beam is being accumulated i. e. unpolarized particles are continuously being fed into the 
beam, 3) the particle input rate is equal to the rate at which particles are being lost due 
to scattering beyond ring acceptance angle, thus the beam intensity remains constant, 
4) increasing the initial polarization of a stored beam by spin filtering, 5) the input of 
particles into the beam is stopped after a certain amount of time, but the spin filtering 
continues. The rate of depolarization of a stored polarized beam on passing through an 
electron cooler has also been shown to be negligible, because of the low electron areal 
densities in cooler beams. 

The principal advantage of this method of polarization buildup is that it can be ap- 
plied to antiprotons. The investigation carried out in this paper will be of use to the PAX 
collaboration and to future projects that hope to produce a high intensity polarized an- 
tiproton beam, or to utilize spin filtering in a storage ring. This work can also be applied 
to the recent idea of polarizing a beam by channeling through a bent crystal [23] . 
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